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Abstract 

Invoking an initial symmetry between the time t and some extra spatial dimension 

y, we discuss a novel scenario where the dynamical formation of the 4-dim branc and 
its cosmological evolution are induced simultaneously by a common t y symmetry 
breaking mechanism. The local maximum of the underlying scalar potential is mapped 
onto a 'watershed' curve in the (t, y) plane; the direction tangent to this curve is identi- 
fied as the cosmic time, whereas the perpendicular direction serves locally as the extra 
spatial dimension. Special attention is devoted to the so-called slinky configurations, 
whose brane cosmology is characterized by a decaying cosmological constant along the 
watershed curve. Such a slinky solution is first constructed within a simplified case 
where the watershed is constrained by y = 0. The physical requirements (and the prob- 
lematics) for a slinky configuration to generate a realistic model are then discussed in 
a more elaborated framework. Prospects for adding matter and radiation fields, their 
localization and their cosmological implications are briefly discussed. 



1 Introduction 



Extra dimensional theories of particle physics and cosmology have received widespread atten- 
tion for more than a decade now. As the typical scale of new physics in phenomenologically 
plausible extra dimensional scenarios is naturally around a TeV, we hope to be able to 
test them at the Large Hadron Collider (LHC) and other contemporary experiments. The 
original idea dates back to the works of Kaluza and Klein (KK) PQ |2] aiming at a unified 
theory of electromagnetism and gravity. In the KK construction [2] all known matter fields 
propagate in the full extra dimensional spacetime and the 4D low energy effective theory is 
obtained by compactifying the extra dimension on a circle and keeping the zero mode of the 
KK tower. More recently, string theory constructions including D-branes and in particular 
the Horava-Witten solution of supergravity [H H] inspired Arkani-Hamed, Dimopolous and 
Dvali (ADD) to introduce the so called Large Extra Dimensions paradigm in order to solve 
the gauge hierarchy problem [SI El [7j. In the ADD model, only gravity is free to propagate 
in the extra dimensional space, while SM fields are confined to a 3 + 1 dimensional hyper- 
surface, referred to as a brane. The ADD construction allows, in principle, for "large" extra 
dimensions already at a mm scale, corresponding to a very small 5D Planck mass. However, 
recent experiments probing the short distance behavior of Newtonian inverse square law have 
already placed upper bounds of (9(40/im) on the size of the extra dimensions for relevant 
realizations of the ADD idea, see for example [HI El HO] . Shortly after, Randall and Sundrum 
(RS) [11] offered an alternative to compactification in the form of a warped extra dimension 
confined between two branes, referred to as the UV and IR brane. In the RS construction, 
it is again only gravity which propagates in the full 5D space time, while the SM fields are 
confined to the IR brane. The warped geometry implies a varying 5D scale along the extra 
dimension which is of 0{Mpi) at the UV brane and of (9(TeV) at the IR brane. By pushing 
the IR brane to infinity, it has been shown [12] that it is possible to localize gravity to the 
UV brane and that deviations from the Newtonian inverse square law become relevant at 
very small length scales of 0{TeV)~^. 

In the ADD, RS and some previous constructions [THl [T3j the branes are fundamental, 
namely they are treated as infinitely thin delta distribution sources. The 4D effective theory 
is obtained by performing a KK decomposition of brane and bulk fields and studying the 
corresponding interactions. Main implications of these constructions are new KK particles 
with masses of 0{TeV) and modifications of the Newtonian potential, both of which are al- 
ready being challenged by LHC results and short range gravity experiments [H]. Depending 
on the theoretical setup, additional experimental constraints can come from precision mea- 
surements of rare decays [15], electroweak precision observables [I6], the Cosmic Microwave 
Background (CMB) anisotropies [17] and more. 

Subsequently, the original extra dimensional setups were modified and extended in the 
context of both particle physics and cosmology. The cosmological evolution of fundamental 
branes was already studied shortly after the introduction of the RS model [TSl [T21 [201 ED 1221 
|23] , exhibiting a non standard cosmological evolution above an energy scale related to the 
brane tension and the curvature of the bulk. In the last decade, various aspects of brane world 
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cosmology have been extensively studied including inflation, gravitational perturbations and 
late time acceleration, with emphasis on possible signatures in the CMB anisotropy spectrum 
[T7j . More exotic models considered the possibility of colliding branes, with fundamental 
branes f2M or within the D3 — D7 string theory inspired constructions |27j . 

If our universe is indeed a 3-brane, coming from a string theory realization of RS or other 
braneworld setups [281 |29l |30], treating it as an inflnitely thin brane suffices for deriving the 
4D particle physics and cosmology. Effects of finite thickness of the brane can be accounted 
for by averaging procedures developed in [311 [321 E31 El]- 

Alternatively to string theory, the simplest underlying dynamics that generates a brane 
may have a field theoretical origin; the brane would correspond to a topological defect 
(domain wall (DW), vortex, etc.) in higher dimensions. The idea that our universe is a DW 
brane in higher dimensional space time dates back to the work of Rubakov and Shaposhnikov 
DW configurations can be generated by a scalar field supported by a double well or 
periodic potential and are able to provide a smooth realization of the brane and the RS 
warped metric [371 EHl EHl SOl IH]. In these setups the width of the DW in the extra 
dimension plays an important role, in particular in the localization mechanisms of gravity, 
gauge and matter fields to the brane [121 HHl HI] . 

What is common to all of these constructions is the maximally 4 symmetric geometry 
associated with the DW brane, allowing for a Minkowski (M4), (iS'4 and AdS^ brane solutions. 
The dS^ case already introduces difficulties in the form of (naked) curvature singularities 
at a finite distance from the DW brane [391 HO], which can be overcome by compensating 
divergences in the energy momentum tensor of the DW scalar [_40j. A dS^ DW solution where 
the singularities are relaxed to be horizons can be obtained numerically [39l B5] ; however, it 
still introduces problems for the localization of fermions and possibly other fields. 

Differently from the maximally 4 symmetric cases, there has been so far only one attempt 
at finding maximally 3 symmetric DW configurations, the latter being relevant for cosmology 
|46j . The resulting solutions are "time shifted" versions of the maximally 4 symmetric 
solutions of pij, yielding a bouncing cosmological evolution opposite to the desired one: 
H[r > 0) < and H{t < 0) > 0, where r denotes a conformal time coordinate, related to 
the cosmological time t by a{T)dT = t [45]. 

The effective cosmology of fermionic and scalar fields in the vicinity of an unspecified 
time dependent DW configuration in 5D was studied recently in [17], as a sequel to previous 
attempts to obtain the Standard Model (SM) of particle physics as the low energy effective 
theory on a Minkowski DW brane in 5D [IH]. In this setting electroweak symmetry breaking 
(EWSB) is driven by the DW scalar itself. A subtle interplay between the bulk scalar fields, 
the DW scalar and gauge fields is needed in order for the SM particle content (Higgs and 
fermions) and the SM symmetry group to be confined to the brane. This construction was 
later extended to a larger Grand Unified Theory (GUT) symmetry group in [49j. For a 
thorough overview of the various aspects of DW configurations (Thick Brane solutions) the 
interested reader is referred to [SO] . 

^Higher dimensional topological defects were also considered around the same time [36], but are not the 
focus of the present work. 
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In this paper we aim at obtaining cosmologically plausible time dependent DW configu- 
rations, following a rather different approach from the above attempts. This tells us that we 
have to adopt a Friedmann-Lemaitre- Robertson- Walker (FLRW), maximally 3 symmetric 
ansatz for the 5D geometry, with crucial implications as will be explained below. 

The paper is organized as follows. In Sec. |2]we elaborate on our approach to the problem 
and define the notion of "slinky" DW configurations. In Sec. [3] we obtain the simplest 
solution with most of the desirable characteristics described in Sec. |2} We then analyze the 
associated scalar field energy densities, and attempt to distinguish between brane localized 
and bulk energy densities; this will provide a better understanding of the effective dark 
energy density induced on the brane. In Sec. |4] we perform a matching of the solution of 
Sec. |3] to an instantaneous RS-like action, to better understand the localization properties 
of the instantaneous bulk and brane cosmological constants. Subsequently, we analyze the 
resulting brane cosmology for early, late and intermediate times. In Sec. |5] we discuss the 
problematic nature of the solution of Sec. |3| due to the non covariance of the resulting scalar 
potential and offer simple remedies and alternatives. Finally, we conclude in Sec. [6] 

2 General Strategy 

As mentioned in the introduction, previous works on extra dimensional codimension one do- 
main wall branes (or thick branes) have so far considered maximally 4-symmetric geometries 
(M4, AdS^, and (iS'4). Consequently, the solutions obtained in this context have limited rele- 
vance for cosmology, which is characterized by a maximally 3 symmetric geometry described 
by a FLRW metric. 

The difficulty in obtaining solutions to the cosmological (FLRW) case stems from the fact 
that all metric coefficients are {y,t) dependent, where y is the extra dimensional coordinate 
and t is the time coordinate. The resulting Einstein and Klein-Gordon equations for the 
evolving scalar field(s), translate into a system of non linear partial differential equations, 
which is extremely difficult to solve even in the most simplified cases, as we shall see below. 

The maximally 4 symmetric DW configurations mentioned above [23 SOI SI] , are charac- 
terized by a scalar potential with two degenerate minima and a (static) kink mode interpo- 
lating between the two minima. The nature of the brane (M4, AdS^ or dS^) is studied in the 
infinitely thin brane limit [IQ], where the DW is squeezed to its central point {y = 0). The 
brane type is then determined from the brane induced cosmological constant and the bulk 
cosmological constant, or more precisely by whether they deviate from the RS fine tuning 
relation 

«:(A;,)V6 + As = 0, (1) 

where A;, denotes the brane induced cosmological constant coming from the kinetic and 
potential terms and A5 is the (negative) bulk cosmological constant, which is identified with 
the minima of the scalar potential. 

What we wish to emphasize is that in order to study the cosmological evolution of a 5D 
time dependent DW configuration, we actually have to study the evolution of the spatial (3D) 
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Figure 1: Schematic view of the watershed in y — t coordinates. The curve 7(7/(5), represents 
the points, which are mapped to the maxima of the scalar potential. The surfaces S|t-=t-o denote 
constant cosmological time slices, where in the most general case the cosmological time is identified 
with the tangent to the curve 7(7/(5), t(s)), and does not necessarily coincide with the 5D time 
coordinate t. The cosmological evolution of such a configuration will be studied by matching the 
instantaneous bulk and brane cosmological constants on each of the S|t-=to hypersurfaces. 



geometry along the y and t directions, simultaneously. This means that in the most general 
case, there is no a-priori reason to consider y = (or any other constant y hypersurface) as 
the position of the brane in the extra dimension. In general, we expect there will be some 
curve in the y — t plane which is mapped by the DW scalar, 0, to the center of the DW 
configuration (0 = 0), which is in turn mapped by V{(j)) to its maximum. Recall that the 
existence of a kink mode implies a Z2 symmetry of the potential which corresponds to two 
degenerate minima |40j . 

The above situation is depicted in Fig. [T| where the curve 7(7/(5), t(s)) represents the 
points in the y — t space, which are mapped to the maxima of the scalar potential, thereby ac- 
quiring the name "watershed" . In the most general case the tangent to the curve 7(7/(5), t(s)) 
serves as the cosmological time, while the normal direction serves (locally) as the extra spa- 
tial coordinate. Namely, at each constant cosmological time slice, we can study the local- 
ization properties of the kink energy density and distinguish between brane localized and 
bulk energy densities. This provides the effective instantaneous brane cosmological constant, 
Al"^^' = + A5. In this way, we actually study the time evolution of a kink induced 

dark energy density on a 3-brane located at (and evolving along) 7(1/(5), t(s)); 7(7/(5), ^(5)) 
is generated by the kink itself. A fully realistic cosmological evolution will necessitate the 
presence of extra sources, for which the most simple example is an (a-priori) brane localized 
perfect fluid with equation of state p = cufp. These additional sources will correspond to 
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radiation and matter densities on the brane, while the kink will contribute the dark en- 
ergy part. Finally, we wish to emphasize that the above prescription is physically plausible, 
only as long as the tangent to the curve 7(1/(5), is timelike for any value of the affine 
parameter s. 

The 5D DW configurations thus generated will be cosmologically plausible if they can 
account simultaneously for a (finite) early time infiationary period and late time acceleration. 
This can be achieved by an induced brane dark energy density, which is extremely large at 
very early times and extremely small at late times. For a better intuition of their time 
evolution, we employ a simple pictorial analogy with a "slinky" spring: " At t = all the 
links of a slinky spring, which represents the DW scalar, are sitting on the maxima of the 
scalar potential between two degenerate minima. As time goes by, the links gradually fall 
towards the minima on the right and on the left, so that at late times, only one link of the 
slinky is interpolating between the two piles of links sitting at the two minima". 

The early time accumulation of links at the maxima corresponds to a huge initial dark 
energy density, decreasing with time at a yet unspecified pace. This energy density drives an 
infiationary period on the DW brane, which terminates when most of the links of the slinky 
hit the minima, or equivalently when the DW kink configuration becomes very thin. The 
late time interpolating link(s) will thus correspond to a small remnant dark energy density, 
which can drive the observed late time acceleration of our (DW) universe. 

Given these requisites, the simplest slinky solution we can think of are scalar profiles, 
nearly fiat along the extra dimension at early times and evolving to step functions at late 
times. The location of the jump or equivalently, the center of the DW configuration, will 
correspond to the position of our brane universe in the extra dimension. Thus, the thin 
brane limit of slinky configurations is achieved dynamically at late times, instead of being a 
mathematical limit of the parameters entering a static DW solution, as in the maximally 4 
symmetric cases [40j. 

A feature that cosmologically plausible slinky configurations should also possess is the 
ability to generate a "brane in time" in addition to a "brane in the extra dimension" . Namely, 
in order for such configurations to account for the creation of the evolving brane universe 
described above, the hypersurface S|r=0! corresponding to zero cosmological time, should be 
distinguished from all other constant r surfaces, and identified with the thin brane limit in 
the (cosmological) time direction. 

A configuration that possesses the above highly nontrivial property will naturally corre- 
spond to the Hartle-Hawking no boundary proposal ^T\ . 

3 Setup and the simplest slinky configuration 

The only recorded attempt at finding time dependent DW solutions is in [36|, where the 
simplest ansatz for the geometry is employed, with a single warp factor depending on time 
and the extra dimensional coordinate. As a first step towards obtaining DW slinky configu- 
rations, we seek for a solution with the desired properties and for the same metric ansatz. 
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We start by defining the setup. The 5D action is given by: 

R 1 



S = / d^xVG 



2k 2 



with metric ansatz 



ds"^ = GAsdx^dx^ = a^{u, r) [dr^ — dx"^] — b'^{uj, T)duP' 



(2) 



(3) 



where k = SvrGs = 87r/M| and G5 is the 5D gravitational constant. The Einstein equations 
for the action in Eq. ^ can be written as follows: 



Rab = 



d 



'A(P0B(, 



-GabV 



(4) 



We explicitly write the components of the Einstein equations in a form analogous to [46j . 
For the (00) component we get: 



62 



[h' + 4h^ -hf]- F + 3H + F{F - H) 



= n 


02 - Iva^ 




3 



(5) 



where ()' denotes differentiation with respect to u and () denotes differentiation with respect 
to r. In addition we have defined H{u,t) = d/a, h{u,T) = a' /a, F{uj,t) = b/b and 
f{uj,T) = b'/b. For the (ij) = 1,2,3) component we obtain (we omit 6ij): 



62 

The (55) component yields: 
62 



H + 2H^ + HF 



(6) 



F + F^ + 2HF 



Ah' + Ah{f -h) = K 



b'^ + -Vb"" 
3 



(7) 



Finally, the (05) component provides a momentum constraint (on 0) of the following form: 

fi-AH' + 3Fh = . (S; 
The (00) and (ij) equations can be combined to yield: 

■n 1 



K 



2H^ + 2HF ~F^ -F-2H 



(9) 



Similarly, if we combine the (ij) and (55) components we obtain 

3 „ 1 62 



0'2 =_(/,/_ h') + 



K a 



F - H) +{F + 2H){F - H) 



(10) 
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Turning to the Klein-Gordon equation (G^^V^Vb0 + dV ld<^ = 0) and using the metric in 
Eq. ^ we get: 

■ dV 

- + {2H + F)<j>- -{Ah - f)<j>' + = . (11) 

Having stated the relevant equations we look for the most simple realization of the slinky type 
configuration introduced in Sec. 2. The first solution we are going to inspect corresponds to 
the simplest watershed possible, a line of constant u that we choose to be a; = without loss 
of generality. Inspired by the time shifted solutions of [IH], we propose a time "proportional" 
kink solution, with a non factorizable dependence on co and r. A simple example of such a 
kink configuration is: 

oc arctan [Awr]. (12) 

It satisfies (j){uj, r = 0) = 0, which we interpret as the initial state of the slinky configuration, 
where all the links of the slinky (0 values along w at a constant r slice) are sitting on the 
maxima of the yet to be determined potential. The time evolution of the above configuration 
dictates that as time goes by more and more links of the slinky will fall towards the minima 
of the potential as approaches its asymptotic value. Although the potential and warp 
factors are yet unspecified, it is clear that a kink solution will require a Z2 symmetry of the 
potential, which corresponds to two degenerate minima symmetric around = 0. For the 



kink configuration of Eq. ( 12 ) the Z2 symmetry translates into refiection symmetry around 
u = 0, under which the kink has to be odd. In addition, since we are interested in the 
localization of gravity to the dynamically generated domain wall brane universe, we should 
consistently seek for warp factors that are peaked at u = and are Z2 even. 

We now consider the time "proportional" generalization of the the maximally 4 symmetric 
warp factors in the spirit of the time shifted solutions of [36]: 

a{uj,T) = — ^ b{uj,T) = ea{u,T) . (13) 

V 1 + A'^U^T^ 

Substituting the above warp factors in Eq. ^ and integrating, we obtain a time propor- 
tional solution for 0: 

arctan [Xur] + /i(a;), (14) 



with integration constant fi{cu). Substituting again Eq. (13) in Eq. (10), we similarly obtain: 

arctan [Xujt] + /2(t), (15) 

with integration constant /2(t). Substituting the solution of Eq. (10) in Eq. ^ and vice 
versa, we realize that fi{y) = f2{j) = const. Setting fi{co) = /2(t) = 0, we substitute the 
above solutions for 0, a and b in Eq. ([6]), solve for the potential V^(0) and obtain: 

T..^N 3A2(r2-eV) . , . 2 ^^ 3(AV - tan^ ^/r^) , ^ . . , 
\/($) = ^^2 ^(l-5sm^$) = ^ — ^(l-5sm2$) 

3(tan2$/u;2_^2^2^2) 
= ^ -^2 ^(l-5sm2<|.), (16) 
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where $ = a/k/S^. The above solution for V{^) is consistent with the rest of the Einstein 
equations in which V appears and the Klein-Gordon equation. We postpone to Sec. 5 the 
detailed discussion of the problems related to an explicit coordinate dependence in the above 
potentiaQ We first elaborate on what it is exactly that we are going to interpret as brane 
cosmology in the above setup. 

Observing the time evolution of (f){uj, r), we realize that when r — t- oo we are naturally in 
a thin brane limit, where the brane is located at w = 0. Inspired by the thin brane limit of 
maximally 4 symmetric domain wall configurations with non vanishing brane cosmological 
constant, A;, [ID], we recall the way in which this cosmological constant is matched. To be 
more specific the brane cosmological constant is due to the energy density of the bulk soliton 
when it is squeezed to the very same brane it generates at a; = 0, once the thin limit is taken. 

3.1 Obtaining the brane cosmological constant in the thick RS 

case 

As an instructive exercise, we review the results of |40j for a maximally 4 symmetric domain 
wall configuration which admits the RS model as its thin limit. The metric ansatz is: 

ds^ = -dy^ + e^f^y) {dt^ - dx""] . (17) 

Considering again a minimally coupled scalar field in the geometry described by the above 
equation, one writes the Einstein and Klein-Gordon equations, which are significantly simpler 
than the time dependent case and obtains the following solution: 

ip = Aarctan[tanh(/3y/2)] V{ip) = A^/S^/S - {A'^l3y8){l + kAVs) sm\2ip/A), (18) 

where A and /3 are constants. The geometry that corresponds to the above kink configuration 
and potential is dictated by e^^y\ where: 

e/(y) = D[cosh{/3y)]~^^'/^^. (19) 

The thin limit of the above configuration is a delicate one in which /3 — )■ oo, A — )■ and 
A^IS is fixed. In this limit / ^ -KA^l3\y\/12, such that /' -> -{KA^l3/12){e{y)-e{-y)) and 
/" — )■ —{KA'^f3/6)6{y). The resulting bulk geometry is that of a slice of AdS^ divided by the 
brane (infinitely thin domain wall centered at ?/ = 0) into two disconnected regions, across 
which there is a jump of the extrinsic curvature. On identifying the degenerate minima of 
the potential in Eq. ( 18 ) with the 5D cosmological constant associated to the thin-limit AdS^ 



space, one obtains: 

V{±7iA/A) = A5 = -KA^py2A f" -2{-K5K/Qf'^5{y). (20) 

The presence of the delta function in /" combined with the (00) component of the Einstein 



equations for the metric in Eq. (17) [10], entail that both ip''^ and V{(p) must contain the 



^Notice that the coordmate dependence renders the potential and thus the whole solution to be non 
a; -s-> T symmetric. 
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same delta function. Consequently, in the thin limit a term e^'^[v9'^/2 + V(y9)] — )■ A^/3/2(5(?/) = 
A.bS{y) is generated in Too- This term corresponds to a brane cosmological constant, which 
automatically satisfies the RS fine tuning relation A5 + = 0. Namely, what had to be 

put by hand in the RS model is an output of the above dynamics. 

An alternative way to derive A;, is by requiring the 5D action / dy^f^C = - J dye^^[ip'^ /2+ 
V{ip)] = (3/«;) / dye^^if" + 2p) to generate the RS action / dye'^f{-Ar^ - AbSiy)) in the 
thin limit. Using this prescription it is clear that A5 comes only from the 2/'^ term and A;, 
comes only from the /" term. The RS fine tuning relation is obviously satisfied by these 
results. 



3.2 Obtaining the brane cosmological constant for slinky configu- 
rations 



In our case we will assume that the time evolution of the kink configuration of Eq. (14), 
also controlled by the dimensionless parameter A, is slow enough that a thin brane limit can 
be taken at every constant time slice, similarly to the maximally 4 symmetric case of [IQ]. 
Then, by matching to an instantaneous RS-like action or inspecting the behavior of Too, we 
determine the bulk {A^{t)) and brane induced (A;,(r)) cosmological constants at each slice 
of constant r. The brane (bulk) cosmological constant is an energy density which is constant 
in the instantaneous 3(4) space dimensions and satisfies p = —p. 

To further elaborate along the lines of the above procedure, we first write the metric 
components along a constant r slice, at say r = tq: 

ds'^l = qf^dx^dx^ = —a^iu, r = Tn)dx^ — e^a^ioj, r = Tn^duj"^ = :z Idx'^ + e^rfw^l, 

\t=to V 5 u/ V 5 Kij 1 + A^CJ^ 

(21) 

where A = tqA and spatial indices j = x, y, z, lo. We note in passing that the metric ansatz 
we chose in Eq. ^ is conformally fiat, but in contrast with the maximally 4 symmetric cases 
in [371 EH ES Sni IH], it does not admit an asymptotic AdS^ limit (with a 5D cosmological 
constant corresponding to the value of the minima of the potential), due to the non-removable 
coordinate dependence in V{(j)). Moreover, although we have written the time and extra 
dimension as r and respectively, we recall that the t — uj symmetric ansatz for the warp 



factors in Eq. (13), is inspired by the solutions to the maximally 4 symmetric case with 
a conformal extra dimensional coordinate |16]. Thus, in order to inspect the cosmological 
time, t, and proper distance coordinate, y, we should in principle rewrite the solution using 
the relations dy = b{uj, T)duj and dt = a{uj, T)dT |^ Provided no singularity occurs in the 
coordinate transformation, we can work with {u, r) coordinates and later transform to {y, t) 
coordinates without loss of generality. 

Considering the instantaneous geometry on each constant time slice described by the 



metric in Eq. (21), we obtain the bulk energy density and match it to a RS-like setup as 



■^We have not specified whether r is indeed a conformal coordinate, but merely stated that the ansatz for 
a and b is inspired by a solution to the maximally 4 symmetric case, where the uj coordinate is conformal. 
This is also true for t in the "time shifted" solutions of HSl. 
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follows. From the action in Eq. or equivalently from the (00) component of the energy 
momentum tensor Too = (1/2)02 ^ (i/2e2)0'2 + a^v{(l)) (where Tmn = {2/^)5C/5G^™ 
and the Einstein Equations take the form Rmn — {^/'2)GmnR = ^^Tmn), we realize that we 
have three separate contributions to the energy density which we label as: f^^/a, ^2^2 and 
fiy(^), where Too = ^</.'2 + ^^^2 + ^v{(j>) and 

O - J_^'2 _ ^r^A^ _ SA^ 

- 2^2*^ - 2k62(1 + A2a;2^2)2 - 2^^2(1 + p^2)2' 

_ 1 • 2 _ 3u^X^ _ 3uj'X' 

" 2*^ " 2k(1 + A2a;2r2)2 " 2^(1 + A2^2)2' ^^3) 

2....^ 3AMr2-a;V)(l-4rW) 3 (a^ - 4AW - A W + 4A We^ 

^^vu) = aV(0) = ^ — K = — ^ . 

2A:e2(l + r2c.2A2)2 2^^2(1 + 3.2,2)^ 

(24) 

In the above equations the second equality is written on E|^=t^o' which we are going to 
characterize each of the above contributions, term by term, to determine the localization 
properties of the associated energy densities and consequently what should be interpreted 
as Ab{To). 

To implement the procedure described in the above paragraph, we perform a transfor- 
mation — 7- y to a proper distance extra dimensional coordinate defined on each constant 



time slice, T,\r=To- Using Eq. (21) we obtain: 



a{u,T = To)du = — zdu = dy =^ ?/ = arcsinh(A(X')/A =^ a; = sinh(A?/)/A. 

VI + A2w2 

(25) 



While the scaling of each term in the energy densities in Eqs. (22)-(24) with the time co- 



ordinate T is determined by the powers of A, their u dependence remains explicit and is 



translated into the dependence on y according to Eq. (25). The thin limit on each constant 
r slice (when tq is simply a finite constant), is taken in principle by letting A — )■ 00. However, 
this limit is not well defined since the A dependence of all contributions is such that most 
of them will naively vanish in this limit. We can in principle overcome this problem by 
modifying 0, V{(f)) and a{u, t = tq) to contain an additional pre-factor A such that the thin 
limit is taken by keeping A^X fixed while letting A and A — ?■ cxd as in [40]. However, such 
a modification of the warp factor, a{u,T = tq), by virtue of Eqs. ^ and (10), will result in 



a non solitonic profile which blows towards infinity and will therefore fails to generate the 
thick brane (kink) we are looking for in the first place. 

On the other hand, by simply observing the u dependence of the same contributions 



in Eqs. (22)-(24), we see that they are all (except for the last term in Eq. (24)) peaked at 



(or around) lo = 0{y = 0). In particular, the common denominator for all contributions. 



D{X,uj) = l/( V 1 + A2(X'2)'^ translates into [sech{Xy)]'^ , which is strongly peaked at y = 
when A — !■ 00. Moreover, the u dependence of f2^/2, Q12 and flv{(j>) is characterized by three 
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distinct behaviours, uo^D and uj'^D, each of which also enters with different powers of 
A. Since u = sinh(A?/)/A, it is clear that the u^D terms will also be peaked around u = 0. 



So, a priori it seems like all of the energy densities in Eqs. (22)-(24) should be interpreted 



as contributions to the brane cosmological constant and do not correspond to any bulk 
cosmo logical constant, except for the fourth term in fly In order to determine the brane 



induced cosmological constant we will integrate the densities of Eqs. (22)-(24) over u on 
S|t-=t-o; to get a result which we are going to interpret as the instantaneous cosmological 
constant on the 3-brane located at a; = 0, at r = tq. Labelling the four terms in Qy of 




24) as n 



1,2,3,4 
V ' 



we realize that n„ 



Qy and 



-A 



2^2 



/(4A^ 



-n'y. 



e_ limits of integration oj — )■ ±oo translate into y — )■ ±oo, as it can be inferred from 

cosh.{Xy)dy we get: 



Eq. (25). Starting from Qm2 and recalling that doj 



dy 



3A2 



2fi;e(cosh®(A?/)) 




16tanh(A?/) 



35A 



JZ-I-— >+oo 



+ 



48A 

35/te' 



(26) 



where the . . . stand for extra terms which are proportional in addition to sech{\yY'^'^ . These 
terms yield negligible contributions when the y — )■ ±oo integration limits are taken. This 
will be especially true when the A — ?■ oo limit will be taken, a subtle case, which is treated 
separately with a different procedure in Sec. |4} Nevertheless, an explicit calculation of the 
integral in this limit will yield the same A (and e) scaling behavior with a different numerical 
pre-factor, which will correspond to the same late time cosmology on the brane. 

2 3 

Turning to ;2 and the related Vty we get: 



du\/~GVt^', 
which implies: 



dy 



3A2e sinh^(Ay) _ ?>X^e / 8tanh(Ay) 
2kA2 cosh^(Ay) ~ 2kA2 



105A 



>+oo 



+ 



duj^G^f 



4A^ 



V 



A2.2 



dujy/GQ 



32A 



8A^e 
35/€A3' 



(27) 



(2^ 



"e" J 35Ke 

Finally, we turn to fiy, which seems to have different localization properties from the other 
energy density contributions. After integration, we obtain: 



duj^Qi 



6eA2 /• sinh^(A2/) , GeA^ / 2tanh(A2/) 



V 



kX^ J cosh (Ay) 



dy 



kA2 



35A 



>+oo 



+ 



24A2e 
35kA3' 



(29) 



Summarizing, we see that while the contributions of Vty and f2^2 cancel each other, the 
identical contributions of VLff,i2 and VLy add up to yield an energy density peaked at ?/ = 0, 
and scales with A. Together with the contribution coming from fiy, this energy density 
is positive and amounts to ~ 64A/(35Ke). As we shall show below the contribution 

coming from f2y corresponds instead to an instantaneous bulk energy density, which is a 
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direct generalization of the bulk cosmological constant, A5, of the maximally 4 symmetric 
case discussed in |10]. Further insight into the role of these energy densities will be gained 
through the procedure described in the next section. 



4 The cosmological constant from matching to an in- 
stantaneous RS-like action 

To get a better understanding of the localization properties of and Q^2^^,2, in particular 

in their natural (yet problematic) late time thin limit, A — )■ 00, we first conveniently rewrite 
the equations of motion in terms of fa{uj, t), defined as: 



fa{y,r) 



In cosh(A?/), 



(30) 

where the last equality is written on a constant r slice according to the coordinate trans- 
formation of Eq. (25). Since 6(0;, r) = ea(a;,r), we have H{u,t) = F{uj,t) = fa and 
h{uj, t) = f{co, r) = Thus, it is straightforward to rewrite Eqs. (|6]), ^ and ( |lO) ) in terms 
of fa{uj, t) and its derivatives and subsequently use these equations to rewrite the action of 
Eq. ^ in terms of fa- The resulting action acquires the following form: 



S 



2e2 



f" f'2 

J a e _i_ J a _p2 

9 J a ' 9 J a 



(31) 

The matching to an instantaneous RS-like action should be performed in the thin (A — )■ 00) 
limit of the domain wall configuration obtained in Sec. |2| On constant r surfaces S 1^=^-0 the 
RS-like action will consist of a bulk cosmological constant term and a 3-brane cosmological 
constant, where the brane is located at ?/ = (or a; = 0). To match we first transform the 
action of Eq. (31) to y coordinates and require that: 

S = 



dr) j d^xduVCe-^f^^'^^^^^ 



f" 

J a 



-fa + 



J a 



12 



9 ^ ci I 9 



f 

J a 



[j dr) j d^'xdyVGe'f-^y^^^^ (^J + 2^^ 



/2 



A— >oo 



( / e^^rfr) 



dyd^X^/g^{-Abulk) + / d^X^/^{-A3-brane) 



(32) 



Using cosh(A?/) )■ e'*'l^l/2 we obtain all relevant quantities on T,\t-=to, in the thin (A — )• 00 

A— >oo 

limit. 

^ a{y) — > 2e~^l^l 



faiy) _^ -My\ 



faiy) 



A— >oo 

dfajy) 

dy 



A— ^-oo 



-Atanh(Ai/) ^ -~X{e{y) - e{-y)), 



A— >oo 



(33) 
(34) 



12 



d'fajy) 



-Xhech\~Xy) -2X6{y). 



A— >oo 



(35) 



We immediately realize that the A — > oo limit of fa{y) in Eq. (32), will generate a brane lo- 



calized energy density proportional to A, ^i°_i,rane ~ 6-^/ ('te); which agrees with the previous 



results obtained by analyzing the various terms in Too in Sec. 3.2 In particular the brane 



cosmo logical constant, A 



f" 

J a 

3— brane'' 



is related to the contributions of the fi^ 



and Qy densities 



of Eqs. (l22j) and by: 



f" 
A— >oo 



Tn, 



^y-5{y)kthrane 



(36) 



The contribution coming from fa{y) will instead correspond to the aforementioned gen- 
eralization of the bulk cosmological constant, A5, appearing in [ID] and Sec, 3.1 These 
contributions are y independent in the A — 00 limit and will be specified later. To obtain 
the explicit form for the fa and fa terms we act with derivatives with respect to r on /a(w, r). 



since the coordinate transformation of Eq. (25) is defined only on constant r surfaces. Once 
the derivatives have been obtained we express them in terms of y on S|t-=t-o which results in: 



fa — 



dfa{uJ,T) 

dr 



I 2|r — T 



A sinh^(A?/) 



1 + A2u;2 



A 



X 



(37) 



Acosh^(A?/) A^oo 

This contribution also corresponds to a y independent (totally delocalized in y) energy 
density in the A — )■ 00 limit. Finally, we inspect the A — 00 limit of fa in an analogous way 
to fa and obtain: 



fa 



d^fa{uj,T) 



2X^X^tu^ 



dr^ 



|S|t = Tq 



^ ^ {my) 



(1 + A2w2)2 

d{-y)r-iOiy)-di-y)r) 



1 + A2w2 

A2 

A2 



A2 
A2 



2tanh^(Ay) -tanh^CXy) 



(38) 



From Eqs. (32), (37) and (38]) we realize that the tanh(At/)^ terms in fa and fa are added 
to each other, leaving us with both tanh^(A|/) and tanh'^(A?/) terms, which will together 
contribute to a y dependent bulk energy density. Their contribution enters with a con- 
formal prefactor e'"^^"-, implying that the associated bulk energy density, Qf2 i- behaves as 
cos\^{Xy) tanh(A?/)'^2>4)_ However, since the metric is strongly peaked at ?/ = in the A — )■ 00 



limit and we care about the contribution of Vt 



to the induced instantaneous brane cos- 



mological constant, the e 2/<' conformal prefactor does not modify the step function limits of 



Eqs. (37) and (38) relevant for our purpose. Notice also that if we integrate over y, the total 



contribution of these energy densities remains finite. The associated bulk contribution to the 
brane induced cosmological constant amounts to A^^;{° = {<o/k)X^/X'^ in the A — )■ 00 limit. 



This result corresponds to the contribution of Vty = (6/k)(A2/A2) tanh (Xy) in Eq. (|24[). 



density (on T,r=To), ^, 



f 

J a 

bulk 



Finally, from Eqs. (32) and (34) the contribution of the f'^ term to the bulk energy 



-(6//t)A2/e2 corresponds to the contribution of fly, which is 



thus interpreted bulk energy term. 
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4.1 Simplest Cosmological scenarios 

The above analysis tells us that in the thin (A — oo) limit the slinky configuration introduced 
in Sec. [2] is equivalent to a brane in a ArfS's-lik^ bulk characterized by the warp factor 
(^{y) = e~^'^'. We have found that the bulk and brane cosmological constants are given by: 

From the above equation, we immediately realize that the RS fine tuning relation is asymp- 
totically satisfied on every constant r slice, in the limit A — )• 00 (;large times) thus supple- 
menting us with a static brane. This happens for any value of the parameter e. Also notice 
that, unlike the solution in [ISj, e = 1 does not correspond to a free solution. At any finite 
A, we have an induced brane dark energy density given by: 

^^^nd^ (Af)' ,sl 6AVI 1\ 6A2 6A2 ^^^^ 

Notice that the brane cosmological constant, A^'' enters quadratically in the RS fine tuning 
relation. This is in exact correspondence to the well known analysis of fundamental brane 
cosmology by Deffayet et al. [TS], [22] , as follows. In the general time-dependent fundamental 
brane case, one solves the 5-dimensional Einstein equations by imposing the Israel matching 
conditions relating the discontinuities in the derivatives of the metric components across 
?/ = 0, to delta distribution sources (see [H]). It turns out that the behaviours of the brane 
sources and the metric components evaluated at ?/ = are independent of the metric solutions 
in the bulk, and obey a modified Friedmann equation where the brane energy density enters 
quadratically. 

Thus, the late time (large A) cosmology of the simplest slinky configuration reported 
in Sec. [2| in the absence of additional sources, is simply driven by a dark energy density 
{^^^'^') siinkyi which wiU consist only of the bulk contributions generated by the fl and fa 



terms in the action of Eq. (32), which decay like 1/t for large values of r. 

To conclude the analysis of the late time (r >> 1) cosmology on the brane we write 
(and solve) the effective Friedmann equation for the evolution of the scale factor of the 
instantaneous 3-brane: 

H,{rf = kA^"- = a,(r) = ""^-^^^ = C,r^\ (41) 

where Hb{T) is the brane Hubble constant and Cb is an integration constant. Thus, the 
late time brane is accelerating with a power law for the scale factor, corresponding to a 
deceleration parameter, q = —abab/o-l = — 1 — —0.6, which is rather plausible if we 

simply treat r as the cosmological tim^ 

'^AdS^-like means that constant r slices look like a slice of AdS^ with different conformal factor. The 
latter smoothly varies with time. 

^In the ACDM model, the observed dark matter and dark energy density parameters are ^2] fijvf — 0.3 
and fl\ ~ 0.7, respectively. Together, they yield a density parameter, qACDM — ^^a//2 — f^A — —0.55. 
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The addition of (a-priori) brane localized sources (matter, radiation), will obviously en- 
able various modifications of the brane cosmological evolution, all of which will be analogous 
to the analysis of [IHl 122] • Needless to say, the latter possibility is highly implausible due to 
the large element of arbitrariness or the lack of underlying dynamics. 

Finally, we comment on the early time cosmological interpretation of the solution reported 
in Sec. |2 Since at early times (A — )■ 0) the kink configuration and warp factor are nearly 
flat (</) — 7- and a(u;, r) — 1), we should expect that the average of the large early time bulk 
energy densities, flj-^ j^, will dictate the cosmological evolution of the 4D geometry obtained 
by integrating out u (or y). 

The main problem with this early time picture is the absence of a brane, or a speciflc 
region in the extra dimension, corresponding to the 4D DW brane universe we wish to obtain 
in the flrst place. It is at this point that we take advantage of the fact that the solution 
reported in Sec. [2] is symmetric in u and r. By performing the same analysis of Sec. |4] on 
constant u slices, we realize that in the limit w — )• oo the scalar fleld, 0, becomes a 

"brane in time" (r-brane) localized at r = 0. Similarly, the fa term in the action (Eq. (31)), 
will become a delta distribution in time for w — )■ cxd , in analogy with Eq. (35) and will thus 
correspond to a brane localized energy density. 

Given the above, the hypersurface S|t-=o will correspond to the center of a thick "r- 
brane", with a brane localized energy density that diverges as a; — )■ oo - thus providing an 
element of creation, which occurs smoothly for any flnite value of u, in the spirit of the 
Hartle-Hawking no boundary proposal. 

To conclude, we flnd the cosmological interpretation of the slinky solution introduced in 
Sec. [2] and studied above to be plausible and problematic at the same time. It is plausible, 
because an extremely simple 5D setup with a single scalar fleld is able to generate a DW 
brane universe with an early and late time accelerating phases, thus addressing directly the 
dark energy paradigm. 

The main problem with the same interpretation stems from the subtle identiflcation of 
early time dynamics and the interpolation between the late time and early time regimes. 
The most natural explanation, which seems to suggest itself, is that our DW brane universe 
is localized towards a; — )■ cxd as r — >■ (which can also be perceived as its creation at this 
point in the extra dimension) and evolves to an inflnitely thin brane at a; = for r — )■ oo, 
which corresponds to a different (non trivial) watershed than the one we naively started 
with, namely the line u = 0. Further study of these problems and the search for alternative 
solutions will be the subject of future publications. 



5 l/(<l>)-Implications and alternatives 

We recall that the slinky solution just obtained is associated with a potential that contains 
an explicit coordinate dependence, thus violating 5D general covariance. As a result, the 
momentum constraint given by the (05) component of the Einstein equation (Eq. (|8|) is not 
satisfled, as it corresponds to 5D conservation of energy and momentum, which is itself a 
result of 5D covariance. 
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Specifically, for the kink in Eq. (14) and the warp factors in Eq. (13), the two sides of 



Eq. |8] for the (05) component take the following form: 

; .rr, ^, SAW 3AW , • , 3AW , 

''-''^+^^''= (l + AVry + (l + AW) * = (TTAWji ■ 

The apparent physical meaning of the violation of the above equation is that the energy 
flux along the 5th dimension (or equivalently the 5 dimensional momentum density) encoded 
in Tffg is not enough to account for the analogous flux associated with Gqs! namely, the kink 



solution in Eq. (14) does not inject enough energy in the 5th dimension in order to support 



the geometry described by the warp factors in Eq. (13). Attempting to identify the extra 



term coming from Gqs with an additional scalar source, results in a new set of equations for 



which the slinky profile in Eq. (14) is not a solution. Importantly, notice that as a; — )• 0, oo 
and/or r — )• 0, oo the violation of the (05) equation vanishes, precisely in those regions 
where our configuration mimics a brane in the extra dimension or in time. In addition, 
the violation of the (05) equation is odd in uj and r and will vanish if we integrate one (or 
both) of them out. It is clear that as long as one considers a potential V{(f)) with an explicit 
coordinate dependence, local properties cannot be inferred from 5D covariance away from 
the asymptotic regions. 

Having stated the problems associated with ^($), we try to better understand the dy- 
namics of the slinky links (portions of the slinky configuration itself) on the terrain of the 
scalar potential. 

From the analysis of Sec. Ill on the S|t-=t-o hypersurfaces, we realize that we can rewrite 



y($) of Eq. (16) covariantly as follows: 



o\2 / \2 2 

y($) = ^ ( ^ - J^(tan2<^>) 1 (1 - 5sin^<|.), (43) 

where $ = -^//t/S^. The above potential is still of a double well form for small values of A, yet 
the location of the minima is different on each constant r slice and overlaps the asymptotic 
value $ = ±7r/2 only in the r — )■ oo limit. In this limit, V{^) evolves to a "chopped double- 
well" shape, where the edges become the location of both the minima and the infinitely high 
walls of V{^), corresponding to the asymptotic values $ = ±7t/2. We depict the form of 
the potential for early, intermediate and late times in Fig. |2} The location of the majority 
of the links will always correspond to the potential in the vicinity of the asymptotic values, 
$ = ±7r/2, except for the r = hypersurface, where all the links sit on the local maxima of 
r($) at $ = 0. 

We conclude that most of the links of the slinky experience a fall from the two walls 
of V{^) (around $ = ±7r/2) towards the minima, the latter approaching the $ = ±7r/2 
values as r — )■ oo. Hence, it is only at late times that we can interpret the scalar profile in 



Eq. (14) as a true kink interpolating between the two (cusped) degenerate minima of V{^), 



as speculated in Sec. 1. 
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Figure 2: The scalar potential on the constant r slices, S|t-=t-o (Eq- (43)), for early, interme- 
diate and late times. We set k, = X = 1 and e = 0.9. Notice the increased depth of the two 
degenerate minima in early and late times and the rise of the maxima in late times. 



To conclude this section, we attempt to get rid of the most problematic feature of the 
solution we have found, namely the explicit coordinate dependence in V{^). The immediate 
remedy in this kind of situations is to treat the explicit coordinate dependence in K($) as 
another scalar field (or another source) and solve the resulting new set of Einstein and Klein- 
Gordon equations, so that the extrema conditions for V{^) remain unchanged. It turns out 
that this method does not work for the solution of Sec. [2| since the addition of an extra 
scalar p^(u;,r) = — e^w^ leads to profiles that diverge asymptotically and/or delocalized 
warp factors. Both features are problematic for the purpose of localization of gravity to the 
domain wall. 



By observing Eqs. (16) and (14), we realize that we can also rewrite the potential in a 



way that preserves 5D covariance, by using derivative couplings, namely: 



- 5 sin $ 

cos^ $ 



(44) 



where $ = a//«730. The resulting potential is merely a correction to the kinetic term, which 
in turn changes the equations of motion of Sec. [2] By varying the modified action with 
respect to G^^, the resulting Einstein equations will take the following form. 



Rab = KdA(j)dB(, 
The (ij) component of the above equation becomes: 



1-5 sin $ 



cos^ $ 



H + 2H^ + HF 



0. 



(45) 



(46) 



Together with the assumption b{co, r) = ea{u, r), it implies that a{u, r) = /(r + ew), bringing 
us back to the "time shifted" solutions reported in |16]. These solutions do not exhibit the 
kind of time evolution of the scalar profile that we required for plausible slinky configurations 
in Sec. [2l 
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We add that a consistent solution for the new set of Einstein and Klein Gordon equations, 
which is still of a kink nature, can only be obtained for the free case, namely when the 
potential vanishes. The resulting solution corresponds to the e = 1 case in 



0(w,r) = arctan[A(u; + r)], 1/(0) = 0, b{uj , t) = ea{uj , t) = —====== . (47) 

a/1 + X\uj + tY 

This is a free solution with a non removable time dependence. The cosmological interpreta- 
tion of such a configuration in the spirit of Sec. |2| will consist of a watershed given by the 
null curve w + r = and the same instantaneous energy densities on each constant r slice. 
Thus, the e = 1 case of [16], describes a brane with constant tension, travelling to the left 
at the speed of light in uj and r coordinates. 

Performing again a transformation to a proper distance coordinate, we get: 

, , arcsinh[A(a; + r)] sinh(A?/) 

a{uj,T)duj = dy ^ y = ^-^ oo + t = (48) 

A A 

implying that the metric in {y, r) coordinates is given by: 

ds^ = -b'^ico, T)du^ + a\u, r) Ut^ - dx^] = -dy^ + J- — - Idr^ - dx^] . (49) 

cosh [Xy) 

The above metric describes a static 5D warped space time which is asymptotically AdS^. 
If we would have transformed to a proper time coordinate, t instead of y, by requiring 
a{u, T)dT = dt, we would have instead gotten a metric warped in time, for which the cos- 
mological interpretation is that of an apparent reversed bouncing behavior as discussed in 



^46j . However, such a cosmological interpretation of the solution in Eq. (47) is not in line 



with the point of view adopted in Sec. |2] and thus irrelevant for our purposes. 



6 Conclusions 

In this paper we have proposed a novel type of time dependent DW configurations in a warped 
5D space time. The purpose is to find a dynamical realization for braneworld cosmology. 
After discussing the initial symmetry of the time t and some extra spatial dimension y prior 
to the dynamics generating the DW brane universe, we have introduced the notion of slinky 
configurations. We require that such configurations satisfy the most obvious properties for 
a cosmologically plausible time dependent DW brane setup, namely an element of creation, 
a finite inflationary period and late time acceleration. 

The simplest configuration we were able to find was based on a time "proportional" ansatz for 
the metric warp factors, corresponding to a conformally flat 5D metric. This configuration 
has most of the desired features, however it introduces an explicitly coordinate dependent 
scalar potential. As part of future steps in this direction, it will be important to clarify the 
role and implications of properties such as non factorizability in {t, y) of the solutions and 
how these properties are related to the possible loss of 5D covariance. 
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The cosmological evolution of the proposed DW brane configuration was studied by con- 
structing instantaneous energy densities associated with the DW scalar and its supporting 
potential on constant time (r) slices, where the potential can be written in a covariant form. 
This procedure allowed to identify bulk and brane induced energy densities (cosmological 
constants). The resulting cosmology was shown to include a smooth element of creation 
(in the spirit of the Hartle-Hawking no boundary proposal [SJJ), a naturally terminating 
infiationary period and a late time acceleration phase. 

An underlying dynamics that accounts for more detailed cosmological features would in 
principle require the addition of matter and gauge fields, rendering the problem even more 
difficult to solve. No satisfactory solution has been found up to date. Another, simplifying 
avenue is to study cosmological perturbations on a given slinky configuration. We plan to 
search for new more satisfactory solutions in the context of slinky configurations and follow 
the proposed avenues in order to probe their cosmological implications. 
The long term purpose of this study is to find a cosmologically plausible and fully dynamical 
DW brane setup, which unifies cosmological and particle physics aspects, and can be used 
for particle physics in the spirit of |17j. A configuration able to satisfy simultaneously the 
constraints coming from both cosmological observations and collider physics experiments is 
certainly highly desirable; it is also appealing the possibility that its building blocks require 
nothing more than General Relativity and Quantum Field Theory in 5D. 
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